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High precision metrology systems based on heterodyne interferometry can measure position and attitude of 
objects to accuracies of picometer and nanorad, respectively. A frequently found feature of the general system 
design is the subtraction of a reference phase from the phase of the position interferometer, which suppresses 
low frequency common mode amplitude and phase fluctuations occurring in volatile optical path sections shared 
by both, the position and reference interferometer. Spectral components of the noise at frequencies around or 
higher than the heterodyne frequency, however, are generally transmitted into the measurement band and may 
limit the measurement accuracy. Detailed analytical calculations complemented with Monte Carlo simulations 
show that high frequency noise components may also be entirely suppressed, depending on the relative difference 
of measurement and reference phase, which may be exploited by corresponding design provisions. Whilst these 
results are applicable to any heterodyne interferometer with certain design characteristics, specific calculations 
and related discussions are given for the example of the optical metrology system of the LISA Pathfinder 
mission to space. 
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1. Introduction 

Metrology systems based on heterodyne interferometry 
have numerous applications in different areas of mea- 
surement technology where high precision coupled with 
excellent stability and robust operation are required. For 
example, they have recently been used in precise angle 
measurements and characterization of beam alignment 
[l| , @ , measurement of thermal deformation and expan- 
sion (dilatometry) , in stellar interferometers for 
astrometric measurements [f|, Q;0> and in gravitational 
wave detectors @ , , [H| ■ Heterodyne interferometry 
(dual frequency) offers several distinct advantages over 
homodyne (single frequency) interferometry: The phase 
is read out from an AC rather than a DC interference 
signal, which facilitates and simplifies the detection and 
read-out process on the one hand, as well as avoiding 
detrimental exposure to low frequency noise on the other 
hand. The usage of Acousto-Optic-Modulators (AOMs) 
in heterodyne frequency generation [l2j has spurred a pe- 
riod of continuous improvements in measurement pre- 
cision over the past decades [T^I , i4 [ , [T^] , reaching the 



picometer level for position and the nanorad level for 
angular measurements, respectively 0,[li|,[i3]. 

Accuracy and stability requirements for space-borne 
gravitational wave detectors are exceptionally demand- 
ing so that noise sources and performance limiting fac- 
tors in the respective interferometers must be well char- 
acterized and appropriate mitigating measures taken. 
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Whilst many discussions in this paper apply to hetero- 
dyne interferometers in general, specific numbers and 
design references are given for the example of the Op- 
tical Metrology System (OMS) of the LISA Pathfinder 
(LPF) space mission [l6j. LPF is a technological precur- 
sor to the Laser-Interferometer Space Antenna (LISA) 
mission which aims to detect gravitational waves with 
interferometry [17| and the OMS is the most precise 
metrology systems qualified for space as of today [9] 10]. 
As a fundamental feature of its design, the phase of a 
reference interferometer is subtracted from the phase of 
the actual measurement interferometer so that common 
mode noise fluctuations, which arise from unstable com- 
ponents and optical path sections shared by both inter- 
ferometers, are suppressed to a large extent. However, 
this only holds for low frequency fluctuations, whereas 
high frequency noise directly couples into the measure- 
ment band. Recent experimental tests observed a dis- 
tinct variation of measurement performance depending 
on the relative phase between the two interferometer 
signals. The calculations and simulations performed 
in this paper demonstrate that the rejection of high 
frequency noise depends on the relative interferometer 
phase, which could explain the experimental data. Sim- 
ple analytical expressions are found for amplitude and 
phase noise rejection which are in excellent agreement 
with the results of Monte Carlo simulations. 

2. Digital Heterodyne Interferometry 

In this section we give a short introduction to digital 
heterodyne interferometry and define concepts and pa- 
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Fig. 1. (Color online) Schematic of a high-precision hetero- 
dyne interferometry scheme consisting of a measurement and 
reference interferometer. The laser output is split into beam 
1, denoted by the green line on the left, and beam 2, denoted 
by the red line on the right, at beam splitter Pn- A common 
clock signal is shared between the heterodyne source and the 
phase-meter. 



rameters relevant to the discussions of subsequent sec- 
tions. 

2. A. Interferometer Schematic 

The basic design schematic (greatly simplified) of a high- 
precision optical heterodyne interferometer, similar to 
the one used in the LISA Pathinder (LPF) mission 0, 
is given in FigfT] In order to facilitate the following dis- 
cussion, the values for parameters applying to the opti- 
cal metrology system of LPF are given in brackets as a 
baseline. 

The output of a stable laser source is split into beam 
1 and beam 2 at beam splitter Pq before the 2 beams 
are shifted by the heterodyne frequency /hct (1 kHz) 
relative to another when passing through the respective 
Acousto-Optic- Modulators ( AOMs) . Beams 1 and 2 are 
then coupled through optical fibres onto an ultra-stable 
optical bench where each beam is split again into sep- 
arate components at beam splitters PI and P2, respec- 
tively. The beam components reflected towards the left 
form the " reference interferometer" and are combined to 



interfere on the corresponding diode. Those components 
continuing along a straight line after the beam splitters 
are combined to interfere in the "position interferome- 
ter" , appropriately named with respect to the variable 
phase which it aims to measure (represented by a posi- 
tion actuator). In a practical system this phase change 
could arise due to the expansion of a temperature sen- 
sitive component or the scan of an uneven surface [l3j]. 
Alternatively, it could originate from the longitudinal 
movement of a reflective test-mass as in the basic oper- 
ating scheme of the LPF interferometer [t| [l(J. Due to 
the frequency difference between beam 1 and 2 the inter- 
ferometer signals Si(t) and S% (i) display an oscillation 
at the heterodyne frequency fhet with a different phase 
offset cf> and similar amplitude A{ n for each interferome- 
ter. 



Sl(t) = A in [1 + COS (27T/hett - 0)] 



(1) 



The interferometer signals are detected by photo-diodes 
and processed by a digital phase-meter which samples 
the signals at / samp (50 kHz) and performs a Discrete 
Fourier Transform (DFT) for a single frequency bin cor- 
responding to the heterodyne frequency /hot- To this 
end, the phase-meter and the heterodyne signal source 
mixing into the AOM must share a common clock to 
which they are synchronized. Phase-meter data, consist- 
ing of the real and imaginary part of a complex phase 
vector Ft are output after each DFT period T, corre- 
sponding to an update frequency / up = 1/T (100 Hz). 
The digital processing steps of the phase-meter to ex- 
tract the signal phase are described in Eq|2] 
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Q(Ft) = S{k- At) sin 
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(2) 



where N =500 is the number of points in the DFT 
and the heterodyne frequency corresponds to bin num- 
ber m (fhet — m 27r/T). The vector components are 
then transmitted to a computer for further processing 
(e.g. calculation of phase and position or angular align- 
ment). Note that the analyses and expressions for com- 
mon mode noise rejection derived in the succeeding sec- 
tions are based on the phase read-out scheme described 
by Eq(2] which is found in many digital phase-meters of 
heterodyne interferometers. 

The computer may also implement digital controllers 
for laser frequency, laser power and Optical Path-length 
Difference (OPD) stabilization [lfj, which will not be 
discussed further in this paper. OPD stabilization was 
found to be necessary in the optical metrology system of 
LPF to mitigate the effects of optical sidebands arising 
from the cross-coupling of RF-signals going to the AOMs 
[l8j . The corresponding actuator for the servo- loop is 
depicted right after the AOM at the bottom of Fig[TJ 
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2.B. Noise Rejection Characteristics 

In order to obtain a highly stable position measure- 
ment, the phase tf> n of the reference interferometer is sub- 
tracted from the phase 4> n of the position interferometer, 
both of which are assumed to be affected by the same 
phase and amplitude noise fluctuations n. Through this 
processing step the noise contributions of optical path 
sections shared by both interferometers are removed en- 
tirely from the difference phase cj> n — if> n . The corre- 
sponding path sections in our schematic of FigJT] are 
those between the original beam splitter Pq and the sub- 
sequent beam splitters Pi and Pi for beam 1 and beam 
2, respectively, including the optical fibres and AOMs. 
These components are known to be quite unstable and 
sensitive to environmental changes such as temperature 
fluctuations so that a removal of the associated phase 
and amplitude fluctuations is critical for high perfor- 
mance. Phase fluctuations occurring in the other path 
sections which are specific to each interferometer may 
be neglected as these sections are located on an ultra- 
stable optical bench (made from Zerodur) with mini- 
mal sensitivity to thermal expansion. The bench only 
comprises optical elements made from fused silica which 
are hydroxide-catalysis bonded[l|| to the Zerodur base- 
plate so that a quasi-monolithic structure with superior 
stability and negligible sensitivity to thermal expansion 
is formed. However, such a noise cancellation scheme 
only applies to fluctuations at low frequencies (such as 
those induced by temperature swings) well below the 
phase-meter output frequency / up (100 Hz). The ques- 
tion arises to what extent high frequency amplitude and 
phase noise (as often seen in optical fibres) couples into 
the measurement band and affects measurement perfor- 
mance. This will be discussed in detail for the remainder 
of this paper. 



3. Amplitude Noise 

In the following we shall assume that a sinusoidal input 
signal is affected by amplitude noise with uniform linear 
spectral density nid [V/ \/Hz]. The corresponding fluctu- 
ations for each sampling step represent distinct random 
variables which are uncorrelated with each other and are 
assumed to have the following basic properties: 



imaginary part) of the complex phase vector: 
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where Hk denotes the fluctuations of the amplitude noise 
(an) at a given sampling step k and A[ n is the amplitude 
of the input signal. Applying some trigonometric sim- 
plifications and orthogonality relations we obtain after 
some algebra for the measured signal phase (f) an : 
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It is now useful to rewrite Equation |4] in a way that 
allows to extract the mean phase <fi from the expression: 



arg < e 



i<f> 



C i2irmk/N 



fc=0 
JV-l 



arg 



2n 



k=0 



NA, 



k c i2-!rmk/N-4> 



(5) 



The noise fluctuations appearing in the sum of EqJS] arc 
tiny compared to 1 (rik C 1) so that the argument of 
the complex number in brackets is given by the imagi- 
nary components alone. We therefore find for the output 
phase-fluctuations 
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and for the mean-square fluctuations 
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We shall determine the impact of the input amplitude 
noise on the phase error at the phase-meter output af- 
ter the signal has been digitized and processed. Quan- 
tization noise shall be neglected in the derivations but 
is considered in the accompanying Monte Carlo simula- 
tions (at the level of the sensor noise described at the 
end of section 2). To this end, we make the following 
ansatz for the real part (and an analogous one for the 



Considering that the number of points N in the FFT 
is given by N=/ samp // up we obtain for the phase- 
measurement noise induced by amplitude fluctuations 
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and normalizing by the measurement band (f up /2) we 
obtain for the uniform linear spectral density of the 
phase noise 



LSD(A^ an ) 



"id 



V2 



(9) 
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As an example, in the optical metrology system of the 
LISA Pathfinder mission the uniform linear spectral den- 
sity of the amplitude noise is given (with power stabi- 
lization loops on) by n\d/Ai a = 10~ 6 Hz -1 / 2 so that 
the resulting phase noise is around 1 /irad / \/Hz. This 
is smaller than ADC internal noise which is on the or- 
der of 10 ^rad/vTIz for 12% of dynamic range (To) ? but 
may become appreciable for small signal amplitudes. In 
Fig 12^ the output phase noise according to Eqj8]is plot- 
ted against signal amplitude given as a fraction of the 
Full Dynamic Range (FDR) for input amplitude noise 
of a fixed magnitude. The red solid line denotes the 
predictions of Eq|8] whereas the open triangles give the 
results of Monte Carlo simulations, where the standard 
deviation was found for 2000 simulation runs per ampli- 
tude value. The excellent agreement between the ana- 
lytical formula and the simulations confirms the validity 
of EqU 

3. A. Common Mode Amplitude Noise Rejection 

We shall now investigate to what extent common mode 
amplitude fluctuations can be rejected between two sig- 
nals depending on the relative phase between them and 
revert to Eq|H]for that purpose. Assuming that the sec- 
ond signal has a phase tp with respect to the first signal 
and, without loss of generality, setting the relative phase 
of the first signal with respect to the phase of the phase- 
meter clock signal to zero, we obtain for the difference 
in fluctuations of the output phase 

An an = A0 an - A^an 
JV-l 



noise level 10 x larger 
at 10% of FDR 
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For the mean-square fluctuations of An an we then obtain 
from EqUU] 
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Fig. 2. (Color online) Impact of amplitude noise on output 
phase fluctuations: (a) Output phase noise is plotted against 
input signal amplitude, given as a fraction of the full dy- 
namic range (FDR), according to Eq(8] (red solid line). The 
results of Monte Carlo simulations are given by the open tri- 
angles, (b) Differential phase noise An an is plotted against 
the mean phase difference <f> — tp between signal 1 and signal 
2. RMS values of the amplitude noise at each sampling step 
are 10 _3 rad relative to the input amplitude. The predicted 
|Curve (EqfTJJ is given by the red solid line, the results from 
Monte Carlo simulations are given by the blue triangles. 
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Without loss of generality we can replace tp by <p — ip in 
EqUU Inserting the previous definition of the uniform 
linear spectral density of input noise into this equation 
we then obtain for the RMS fluctuations of the differen- 
tial phase between signal 1 and signal 2 
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Comparing with EqjS] we find that the output phase 
noise of the signal difference is increased by a factor of 
2sin[(0 — if))/ 2] with respect to the output phase fluc- 
tuations of a single signal. This relationship is depicted 
in FigfJb, where the red solid curve is described by the 
analytical formula of EqfTJ] and the blue open triangles 
correspond to the results of Monte Carlo simulations. 
For these simulations the difference angle was varied by 
360 degrees in steps of 6 degrees, whilst for each angle 
the standard deviation of the measured phase differences 
was found from 2000 simulation runs. 
We find excellent agreement between the two methods. 
Depending on the relative signal phase, the noise in the 
phase difference may increase or decrease, reaching a 
minimum for both signals in phase and a maximum when 
the two signals are out of phase by 180 deg, as may in- 
tuitively be expected. When the 2 signals are 90 deg out 
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of phase, the noise level for the difference signal reaches 
a value of y2 times the one for a single signal. Note that 
this is the same noise level as for the subtraction of two 
signals with uncorrelated noise sources. 

4. Phase noise 

Phase noise on the input signal presents the second ma- 
jor source for fluctuations of the output phase. In the 
following we shall derive an analytical expression for the 
output phase fluctuations after two input signals affected 
by phase noise have been processed by the phase-meter 
and their phase difference was found. We make again the 
following ansatz to describe the phase-meter processing 
of one input signal 
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where n k denotes the phase fluctuations of the input 
phase at each sample step k. As in the case of amplitude 
noise, the phase fluctuations are assumed to be uncorre- 
lated with properties described by EqJ3] As phase noise 
is contained within the argument of the input signal, it 
is more difficult to treat mathematically than amplitude 
noise and the derivation is somewhat lengthy so that we 
restrict ourselves to an overview over the major steps. 
The input signal arguments are expanded according to 
basic trigonometric relations and sine and cosine terms 
containing the phase noise (pn) fluctuations are approx- 
imated by the following relations 



cos(0 + n/.) = cos 4> — sin ( 
sin(0 + Wfc) = sin <p + cose 



We then find after some algebra 
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Eg uat ion 1151 describes the phase- meter output phase for 
a signal with given phase noise and is used in the next 
step to calculate the common mode noise rejection be- 
tween two signals. 



4. A. Common Mode Phase Noise Rejection 

We are interested in finding the noise rejection in the 
phase difference between two signals, one of center fre- 
quency 4> and the other of center frequency ip relative 
to the phase-meter reference clock. Both signals are af- 
fected by the same phase fluctuations and we aim to 
determine the noise rejection from the differential phase 
0pn — V'pn ■ The following trigonometric identity is useful 
for this purpose 

tan(</> -ib) = — ^ t (17) 

cot <p + tan ip tan <p + cot ip 

Inserting the definition of Eg J15I for tan(0 pn ) and an 
analogous one for tan(-0 pn ) into Eq[l7l and neglecting 
terms to second order or higher in the noise fluctuations 
rife, we find after some lengthy algebra 

tan(0 pn - ip pn ) = 

sm(ip — ip) 

cos(</> — ip) + 2a cos((p + ip) + {P — 7) sin(0 + ip) 

(18) 

Defining abbreviations for the phase difference x — 
<p — ip and the sum of all the small noise quantities 
e = 2a cos(0 -(- ip) + (/3 — 7) sin(0 + ip) we may rewrite 
EqfTSlas follows: 



f(x,e) = (p pn - -0 pn = arctan 



SIM 



~ x — e sin ir, (19) 

where for the last step we performed a Taylor expansion 
to first order in e of f{x, e). In conclusion, for the mea- 
sured phase fluctuations which are induced by the phase 
noise affecting both signals common mode we find after 
some simplifications: 

An pn = 



"pn 



V'pn) ~ ((p-tp) 



-esm 



= — sin 



N 



N 



sm 



N-l 

ip) cos(0 + ip) E sin I 
fe=o ^ 



N-l 

k=a ^ 



b-1>) 
( 47rmA;\ 



N 



(<p — ip) sva((p + ip) 



ir) nk 



(20) 

It is interesting to note that the fluctuations are only 
fully suppressed if the differential phase between the two 
input signals is exactly zero, i.e. <p — ip — 0, whilst 
the fluctuations are maximal at an amplitude of e when 
the two signals are out of phase by 7r/2. In the next 
step we determine the mean-square value of the noise 
fluctuations and obtain after discarding the cross-term 
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Fig. 3. (Color online) The differential output noise An 
is plotted against the mean phase difference <f> — ip between 
signal f and signal 2: (a) Phase noise only is applied to 
the input signal with RMS values at each sampling step of 
10 _3 rad. The predicted curve (Eq[22} is given by the red 
solid line, the results of Monte Carlo simulations are given 
by the blue triangles, (b) Phase noise ofRMS2xlO" 3 radand 
amplitude noise of RMS 10 -3 are applied. The open triangles 
denote the exact results from Monte Carlo simulations. The 
predicted curves if either of the noise sources is applied alone 
are given by the red and black dotted lines, respectively. 



Defining an input phase noise of uniform linear spectral 
density n\d in units of rad/ V Hz] we obtain from Eq|21l 



(Ar 
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(22) 
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Equation [231 is one major result of this paper. The 
analytical results have been verified by comparison 
against Monte Carlo simulations, which is depicted 
in Figj3ji. Excellent agreement is found between the 
analytical prediction of Eqj22] (red solid line) and 
the output data from the simulation (blue triangles), 
where the standard deviation of 2000 simulation 
runs per angle are plotted. The output phase noise 
increases according to a sine law and reaches a maxi- 
mum at 90 deg phase difference between the two signals. 



5. Comparison to experimental observations 

The situation becomes more complex if both, amplitude 
and phase noise, affect the input signals, which is de- 
picted in Fig(3]D- In this case, the resulting differen- 
tial output noise follows neither the prediction for phase 
noise given by EqJ^U nor the one for amplitude noise 
given by Eq|12] but exhibits certain characteristic fea- 
tures of both. If either source dominates but one does 
not know which one, it is easy to distinguish between 
the two experimentally due to the different extrema and 
periodicity of the output noise as a function of cf> — ip: 
The impact of phase noise is minimal for <fr — ip = ir and 
has a period of 7r, whereas the impact of amplitude noise 
is maximal at <f> — if) — n and has a period of 2n. 
During recent test campaigns of the LPF optical metrol- 
ogy system seemingly random changes of the measure- 
ment noise floor (in certain spectral regions) were ob- 
served in between measurements, which upon closer in- 
spection seemed to correlate with changes in the relative 
phase 4>—tp between measurement and reference interfer- 
ometer (see e.g. [I3|). As the relative phase depends on 
the test-mirror position which is variable, this leads to 
a correlation between the measured noise floor and the 
mirror position. These observations are in agreement 
with the predictions made in this paper for the effect 
of high frequency phase noise but more systematic mea- 
surements would have to be performed to clearly confirm 
such a relationship. It should be pointed out that fre- 
quency mixing and nonlinearities due to imperfect sepa- 
ration of polarization states may lead to per iodic errors 
in polarizing heterodyne interferometers [2Q|. However, 
the optical metrology system of LPF is non-polarizing 
and such errors can therefore be ruled out as a possible 
explanation of the measured effects. 
Interestingly, it has also been pointed out that rela- 
tions Eq(22] and Eq02] are similar to equations previ- 
ously obtained to describe the impact of optical side- 
bands spaced around the center frequency at multiples 
of the heterodyne frequency [l8j]. The detrimental effects 
of those were overcome by inclusion of an OPD actua- 
tor for phase stabilization, as depicted in the system 
schematic of FigJlJ 

6. Conclusion 

We have derived expressions describing the coupling of 
high frequency amplitude and phase noise of an input 
signal into phase fluctuations of the digital phase-meter 
output of a heterodyne interferometer. In a next step we 
calculated the expected common mode noise suppression 
between two input signals and found a dependency on 
sin(</> — -0) for phase noise and on sin[(0 — VO/2] for am- 
plitude noise, where (j) — ip describes the phase difference 
between the two input signals. All analytical derivations 
were complemented by numerical Monte Carlo simula- 
tions which were found to be in excellent agreement and 
to validate any approximations which were made. Re- 
cent experimental observations of correlations between 
measurement phase and noise floor in the optical metrol- 
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ogy system of the LISA Pathfinder mission could possi- 
bly be ascribed to the effects discussed in this paper. 
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